Quasimonoenergetic and low emittance ion bunch generation from ultrathin targets 
by counterpropagating laser pulses of ultrarelativistic intensities 
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A new method for generation of quasimonoenergetic and low emittance fast ion/nuclei bunches of 
solid densities from nanotargets by two counterpropagating laser pulses of ultrarelativistic intensities 
is proposed, based on the threshold phenomenon of particles "reflection" due to induced nonlinear 
Compton scattering. Particularly, a setup is considered which provides generation of ion bunches 
with parameters that are required in hadron therapy. 
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Acceleration of ions with superstrong laser beams has 
attracted broad interest over the last years conditioned 
by a number of important applications, such as inertial 
confinement fusion [l|, isotope production Q, hadron 
therapy 0, and etc. Most of these mechanisms are based 
on the indirect processes, when ions are accelerated by 
the space-charge fields induced in the laser-target inter- 
action process. There are several regimes of ions accel- 
eration depending on the laser intensity: Target Nor- 
mal Sheath Acceleration regime Coulomb Explosion 
regime Radiation Pressure Dominant regime Q, and 
the shock wave acceleration mechanism [7[. Typically 
in Target Normal Sheath Acceleration regime, the accel- 
erating gradient reaches of 10 12 V/m, but due to short 
interaction range ~ 1 /jm it accelerates the ions to tens 
of MeV energies. The energy distribution of accelerated 
ions obtained by the laser plasma interaction is usually 
exponential with almost 100% energy spread up to a 
cutoff energy. Only within the last years became pos- 
sible to produce monoenergetic ion beams from targets 
using special technique In these experiments up to 
60 MeV protons was observed. In Radiation Pressure 
Dominant regime, with circularly polarized laser pulses 
when heating of electrons is reduced, it is expected to 
achieve GeV/nucleon energy range at the ultrarelativis- 
tic laser intensities [§]. 

The use of mentioned indirect processes for ions accel- 
eration, because of the lack of directly accelerating laser 
fields of required intensities, has some undesirable fac- 
tors, such as uncontrollable beam parameters, etc. The 
state-of-the-art laser systems are capable of generating 
electromagnetic pulses with intensities exceeding on sev- 
eral orders the threshold of relativism for electrons. For 
heavy particles, laser intensities, at which, e.g., an ion 
with atomic mass number A and charge number Z be- 
comes relativistic, is defined by the condition S > 1, 
where 2 = ZeEX/ Am u c 2 is relativistic dimensionless pa- 
rameter of a wave-particle interaction (e is the elemen- 
tary charge, c is the light speed in vacuum, and m u is the 
atomic mass unit m u ~ 1.66 x 10~ 24 g ) and represents 
the work of the field with electric strength £ on a wave- 
length A (A = A/27r) in the units of particle rest energy. 



Laser intensities, at which relativistic effects become im- 
portant for an ion (S = 1 -threshold value) can be esti- 
mated as: I r = A 2 Z- 2 x 4.55 x 10 24 W cm~ 2 [A//mf 2 - 
The availability of such intensities now is in the scope 
of ELI project and it is of interest the mechanisms 
for direct laser acceleration of heavy particles. The spec- 
trum of direct acceleration mechanisms of nuclei/ions by 
a single laser pulse is very restricted, since one should 
use the laser beams focused to subwavelength waist radii, 
or use subcycle laser pulses. Direct acceleration of ions 
by petawatt laser beams focused to subwavelength waist 
radii has been explored with radially polarized lasers 
llj . where ions energies cover the domain of applica- 



tion in hadron therapy Q. However, the realization of 
this scheme with radially polarized laser pulses of re- 
quired ultrarelativistic intensities is problematic. Hence, 
it is reasonable (for low energy spreads and emittances 
of accelerated ions bunches) to consider coherent pro- 
cesses of laser-particle interaction. Among those induced 
Compton process is practically more effective. Thus, in 
this process a critical value of intensity exists 12] above 
which the nonlinear threshold phenomenon of particles 
" reflection" occurs from the formed moving wave-barrier 
which opens a principally new way for particles acceler- 
ation on the short distances -even shorter than a wave- 
length Tijj . Note, however, that this phenomenon of crit- 
ical field, in principle, takes place for plane waves. How 
the critical-ficld-cffect will change the interaction dynam- 
ics for strongly nonplane -tightly focused ultrarelativistic 
laser pulses, this is a problem for investigation in an each 
specific case, for concrete pulse-configurations (see, e.g. 
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In this letter we propose a new method for genera- 
tion of solid density, monoenergetic and low emittance 
ions/nuclei bunches of ultrashort durations. Accord- 
ingly, we will study laser acceleration of ions/nuclei in 
the induced Compton process with strongly nonplane 
ultrashort pulses, taking into account the specific fea- 
ture of the critical field for generation of fast ion beams 
with low energy spreads and emittances. For realiza- 
tion of the latter we propose to focus two counterpropa- 
gating laser beams of different frequencies and ultrarel- 
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ativistic intensities onto an nanoscale-solid-plasma tar- 
gets with Z/A = 1; 1/2 (for protons and nuclei, respec- 
tively). In particular, we consider the specific range of 
energies that is of great interest to hadron therapy 
with quasimonoenergetic-monochromatic ion beams. 

At first we consider the classical dynamics of an ion in 
vacuum at the interaction with two counterpropagating 
plane waves of carrier frequencies uj\, uj2 (let uj\ > W2), 
wavenumbers ki = {oJi/c, 0,0}, k2 = {—1^2/0, 0,0}, and 
slowly varying electric field amplitudes E\ (ti), E 2 (t 2 ) 
(ti =t — x/c, T2 = t + x/c). Both waves are assumed to 
be linearly polarized along the OY direction: 

El )2 (x, t) = {0, Ei <2 (n, 2 ) COS Wi, 2 Ti,2, 0} . (1) 

Dynamics of an ion in the resulting electric E = Ei + E2 
and magnetic H =x x (Ei — E 2 ) fields is governed by the 
equations: 

dU - eZ (e\ UxU ) d 1 = _eZ_n 1 E 

dt Am u c \ 7 J ' dt Am u c 7 

Here we have introduced normalized momentum II = 
p/(Am u c) and energy 7 = y/1 + II 2 (Lorentz factor). 
When the particle initial transversal momentum is zero, 
and the waves are turned on/off adiabatically at t — > ^foo, 
from Eq. @ for transversal momentum one can obtain 
(coordinate z is cyclic, so II 2 = const = 0): 

Hy (Si, 32) = Hi (n) sinwi-ri + H 2 (r 2 ) sinu;2T2. (3) 

With the help of Eq. ©, one can obtain equa- 
tions for longitudinal momentum and energy, which in- 
clude four nonlinear interaction terms: two of them 
are proportional to H 2 2 ( r i,2) sin2wi j2 7"i i 2 and describe 
interaction with separate waves that can not provide 
real energy change for ion. The term proportional 
to Hi (t\) H2 (T2) sin (wiTi + W2T2) describes interaction 
with the fast interference wave making no contribution 
to real energy exchange too. This term is responsible 
for particle-antiparticle pair production from Dirac vac- 
uum [IBj . The resonant interaction of the ion for ac- 
celeration is governed by the slowed interference wave 
Hi (ti) H2 (T2) sincD (t — x/cj3 p h), where w = wi — ui2 and 
/3 p h = uj/ (u!i + LU2) < 1 are the frequency and phase 
velocity of the slowed wave. Hence, keeping only this 
resonant term, one can obtain the following integral of 
motion in average: 

7 - fiphJlx - 7o - Pph^ox- (4) 

From Eqs. (UJ), ([3]) and dispersion relation 7 = \/l + II 2 
one can see that for the certain values of Hi. 2 (which 
we call critical) the following relation for an ion average 
transversal momentum in the field may be satisfied: 

(Hi,H 2 ) > 70 \/3 P k - (1 - P 2 ph y 1/2 , 
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FIG. 1: The final scaled kinetic energy versus the initial posi- 
tion of the ion xo in units of wavelength Ai . The arrow shows 
position of the target layer. 

at which the slowed interference wave becomes a poten- 
tial barrier causing ion reflection from this moving wave- 
barrier. After the interaction (Hi ; 2 = 0) for the reflected 
ion final energy we have: 

7/ ~ 70 + 2 lQ (3 ph (f3 ph - (3 0x ) (1 - 0> ph ) - 1 . (5) 

Formula ([5|) shows that ion acceleration depends nei- 
ther on the fields' magnitudes nor the interaction length. 
Namely this feature of reflection phenomenon is used 
here for generation of monoenergetic and low emittance 
ions/nuclei bunches of ultrashort durations. 

Then for actual strongly nonplane and supershort ul- 
trarelativistic laser pulses of certain configurations the 
problem is solved with the help of particle-in-cell (PIC) 
simulations. Here we report on the results of the 2D3V 
PIC simulations of counterpropagating waves interaction 
with nanolayers. We use the code XOOPIC, which is a 
relativistic code based on PIC method [16]. Then we con- 
sider two type of ion targets: ionized hydrogen 1 H + and 
fully ionized carbon 12 C 6+ . The simulation box size is 
40Ai x 20Ai in xy plane. Number of cells are 4000 x 200. 
The total number of macro-particles is about 4 x 10 4 for 
hydrogen target and 2.4 x 10 5 -for carbon target. Tar- 
get is assumed to be fully ionized. This is justified since 
the intensity for fully ionization of carbon is about 10 19 
W/cm 2 , while we use in simulations intensities at least 
on five orders of magnitude larger. So, the target will 
become fully ionized well before the arrival of the pulses' 
peaks. Electrons and ions are assumed to be cold in the 
target, T e =T l= 0. 

The laser pulses have profiles sin 2 (nt/Ti^) with pulse 
durations 7i = 5Ai/c, T2 = 3A2/C and Gaussian trans- 
verse profiles: exp(— r 2 /w\ 2 ), with the waists W\ = 6A1 
and W2 = 3A2. The carrier-envelope-phase of the lasers 
is set to zero, so that the electric fields' maximums are 
at the pulses' centers. The first laser [oj-y) is introduced 
at the left boundary and propagates along the x axis 
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FIG. 2: (Color online) Ion layer blown out from a thin foil 
by counter-propagating waves. The (x — y) phase space for 
protons at instants (1) 15T 2 , (2) 18T 2 and (3) 21T 2 . 




FIG. 4: (Color online) The angular distribution of accelerated 
protons versus the transversal position y. 
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FIG. 3: (Color online) The final kinetic energy distribution 
(at instant 21Tb) of accelerated protons versus the transversal 
position y. 
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FIG. 5: (Color online) Density distribution of protons at in- 
stant 21T 2 . 



from left to right, is focused at the target layer. The 
second laser (^2) is introduced at the right boundary 
and is also focused at the target layer. We consider case 
Ai = A2/2 when A2 = 800 nm which for the phase veloc- 
ity of the slowed wave gives: j3 p h = 1/3 and according to 
Eq. ((©) for reflected particle normalized kinetic energy 
gives: jf — 1 w 0.25. Before PIC simulation we have 
cleared up the role of initial conditions. It is evident 
that for the short laser pulses, depending on the initial 
position of an ion, reflection will take place with various 
longitudinal velocities (acquired in the fields), resulting 
to different final energies. This is well seen from Fig. 1, 
displaying the role of initial conditions: the final energy 
versus the initial position xq of an ion for plane waves, at 
vo = 0. The intensity of effective slowed wave is above 
the critical point: Si = 0.6, S2 = 0.3. Numerical results 
show that this model is almost accurate in predicting the 
final energies of reflected particles. Thus, taking into ac- 
count the result for PIC simulations, the target layer is 
placed at the distance of 15.2Ai from the left boundary, 



and the first laser is introduced with time delay to start 
pulse: t de i = 10Ai/c. 

In Fig. 2 ion layer blown out from a thin 1 H + tar- 
get of thickness 8 nm by counter-propagating waves, is 
shown for various instants. Hear n e — n c , where as 
a measure of density we take critical plasma density 
n c = 1.74 x 10 21 cm~ 3 calculated for 800 nm laser. Note 
that electrons are escaped from the target before the 
arrival of the pulses' peaks. At the 21T 2 the ions are 
already free. The final kinetic energy and angular 
<j> = tan -1 (H x /H y ) distributions of accelerated protons 
versus the transversal position, are shown in Figs. 3 and 
4, respectively. The number density of protons is dis- 
played in Fig. 5. The results for carbon foil of density 
n e — 80n c and 4 nm thickness are shown in Figs. 6, 
7, 8. Estimates for transverse (e t ) and longitudinal (e;) 
omittances of the ions bunch show that for protons et 
< 0.1-7T mm mrad and e; < 10~ 8 eV s, while for carbon 
ions e t < O.ltr mm mrad and e/ < 10~ 7 eV s. These 
results are on several orders of magnitude smaller than 
their counterparts in conventional ion accelerators. The 
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FIG. 6: (Color online) The final kinetic energy distribution 
of carbon ions versus the transversal position y. 
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FIG. 7: (Color online) The angular distribution of accelerated 
carbon ions versus the transversal position y. 



energy spreads are: 8£/£ ~ 10~ 2 . 

In conclusion, owing to the threshold character of par- 
ticle reflection phenomenon, efficient generation of quasi- 
monoenergetic and low emittance fast ion bunches of 
solid densities from nanotargets by counterpropagating 
laser pulses of ultrarelativistic intensities has been put 
forward and demonstrated with PIC simulations. Partic- 
ularly, a setup is considered which provides generation of 
ion bunches in the range of such energies and low spreads 
that may have application in hadron therapy. 
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